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In this article, Einstein-Maxwell space-time has been considered in connection to some of 
the astrophysical solutions as previously obtained by Tolman (1939) and Bayin (1978). 
The effect of inclusion of charge into these solutions has been investigated thoroughly 
and also the nature of fluid pressure and mass density throughout the sphere have been 
discussed. Mass-radius and mass-charge relations have been derived for various cases 
of the charged matter distribution. Two cases are obtained where perfect fluid with 
positive pressures give rise to electromagnetic mass models such that gravitational mass 
is of purely electromagnetic origin. 
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1. Introduction 

The so-called Schwarzschild interior solutions obtained by using Einstein's field 
equations corresponding to various spherically symmetric static perfect fluid distri- 
butions usually suffer from the well-known problem of singularity. An uncharged 
incompressible fluid sphere of mass m cannot be held in equilibrium below certain 
radius a — 9m /A and even demands a larger value for a related to physically reason- 
able equation of state 1 . One way to overcome this singularity due to gravitational 
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collapsing of a spherically symmetric material distribution is to include charge to 
the neutral case. It is observed that gravitational collapse can be averted in the pres- 
ence of charge where gravitational attraction is cou nter balanced by the electrical 
repulsion in addition to the pressure gradient I2C3H 

But even then many questions came up regarding the stability of the charged 
sphere and also about the amount of charge. Bonnor^ worked on this type of charge 
inclusion model and showed that a dust cloud of arbitrarily large mass and small 
radius can remain in equilibrium if it has an electric charge density related to the 
mass density by a = ±p. According to Stettner^, a fluid sphere of uniform density 
with a net surface charge is more stable than without charge. However, Glazer^by 
considering radial pulsations showed that Bonnor ^ model is electrically unstable. 
He also explicitly established the effects of electric charge upon dynamical stabil- 
ity^. The work of Whitman and BurchlS, for arbitrary charge and mass distribution 
showed how charged analogue cases give more stability. T hey also applied the pul- 
sation equations to the char ged solution of Pant and Sah which represents the 
charged analogue of Tolman HI solution of type VI, and got unsatisfactory results 
in connection to the boundary condition which is incompatible with densities and 
pressures, de Felice et al. ^ proposed a model for charged perfect fluid and con- 
cluded that inclusion of charge inhibits the growth of space-time curvatur e w hich 
has a great role to avoid singularities. However, very recently Ray et al. and 
Ghezzi >^ have studied the effect of electric charge on compact stars assuming the 
charge distribution is proportional to the mass density. Both the group have argued 
that with the huge amount of charge and strong electric field these type of stars 
would be unstable and eventually would form charged black holes. 

Therefore, starting from 

the very well-known Weyl-Majumdar-Papapetrou | 14 | 15 | 16 [ type general relativis- 
tic charged dust solutions, via Bonnor ^ and several others (some of which have 
been mentioned in the above introductory discussions), one can arrive at the realm 
of charged compact stellar objects. But a genuine question can be raised in this 
context that: what about the normal stars, do they carry electric charge and if at 
all, what will be the stability of their configurations? As such regarding the effect of 
charge an d possibility of holding charge by the stars are not unavailable in the past 
literature | 17 | 18 | 19 | 20 | , As a continuation of such study Shvartsman^H argued that 
while astrophysical systems are usually thought to be electrically neutral, this may 
not be always true in the real situation. His analysis is based on the exchan ge p ro- 

991 

cesses between stars and the surrounding medium. Very recently, Neslusan has 
reminded about the existence of global electrostatic field of the sun and other nor- 
mal stars. He has given a general charge-mass relation q r — [27reoG(m p — m e )m r ]/q 
where eo is the primitivity of vacuum, G is the gravitational constant, m r is the 
stellar mass in the sphere of radius r, m p and m e are respectively the mass of proton 
and electron having charge q, while q r is the global electrostatic charge inside the 
star. He got it to be q r — 77.043m r when q r in Coulomb and m r in solar masses 
corresponding to an ideally quiet, perfectly spherical, non-rotating star. 
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So, motivated by all the above facts regarding the charged sphere in co nnec tion 
to normal stars, we have considered here charge analogue of Tolman VI anc j 
Bayinl^Sltype solutions which represent very i mportant class of astrophysical solu- 
tions. However, in a series of work Ray and Das | 24 | 25 | 26 | fr ave already been studied 
some specific aspects of these type of solutions in different context. In the present 
work we analyze different physical properties of the static charged stellar model. 
We show that some cases provide electromagnetic mass model (EMMM) where all 
the physical p aramet ers, including the gravitational mass, vanish due to vanishing 
charge density BZEE it has also shown that EMMM do exists even with positive 
pressure which clearly contradicts the observation done by Ivanov ^ that EMMM 
are allways associated with repulsive pressure. The effect of charge inclusion in these 
EMMM, along with other solutions, has been investigated thoroughly in connection 
to mass-radius and mass-charge relations. The nature of fluid pressure and mass 
density throughout the fluid sphere have also been discussed. 



2. The Einstein-Maxwell field equations 

Let us consider a static spherically symmetric matter distribution corresponding to 
the line element 



ds 2 



A 2 dt 2 



B 2 dr 2 - r 2 {d9 2 + sin 2 6dcj) 2 ) (1) 

where A and B are function of the radial coordinate r only. 

Then the set of Einstein-Maxwell field equations, in the co-moving coordinates, 
for the above line element may be explicitly written ^31 as 
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where the total charge within a sphere of radius r, in terms of the 4-current J 4 , can 
be given by 



q = 4ir J°r 2 ABdr. 
Jo 



(5) 



3. The solutions to the field equations 
3.1. Bayin's class of solution 

Now, assuming A'/Ar — C(r) and then equating (3) and (4), we get Bernoulli 
equation for B(r) and C(r) as follows 
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Therefore, choosing q(r) = Kr n we get the solutions forp(r) and p(r) corresponding 
to different values of the parameter n (for detail results vide ref. . 
Case I: For n = 1 

K 2 1 f a + 3air 



(ao + air)r 



, 2r 4r 2 L> , 
C C 2 r 



- 4A ' V (^'"7-^ + G ')}' (7 » 

K 2 4/1 1 \ AK 2 12 , D 

- 8 ^(4 h V-c^ + Gl )- (8) 

where Co and G\ are the constants of integration. Here and in what follows D has 
been substituted for (Co + 2r). 



Case I: For n = 3 
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where G2 is the constant of integration. 
Case II: For n — 1 
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where C\ is the constant of integration. 
Case II: For n — 3 



' [ 2{C 1 r 2 + W 2 -(l-2W 2 K 2 )r i } 1/2 
r 2 + Ci] +3K 2 r 2 , 

8tt P = fi(5r 2 - 3Ci) - llK 2 r 2 . 



Case III: For n = 1 



8tt P = — [K 2 -l + (2K 2 - l)(C 3 r - 2)C 3 r] 
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where C 3 is the constant of integration. 
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3.2. Tolman VI type class of solution 

In their work Ray & Das considered the solution of Pant & Sah t discuss 
about the electromagnetic mass model (EMMM) while in another paper ^3 they 
have discussed the role for equation of state regarding the EMMM. However, in the 
present work we would like to find out some other features of the solutions, specially 
the condition to have positive pressure. The solutions given by Pant & Sah 
corresponding to the metric (Eq. QJ, where e"/ 2 = A, e x l 2 = B and K = E 2 /8n in 
the present mathematical nomenclature, become 

P=7i[l-c(n-l) 2 ], (17) 
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The above set of solutions, in view of c, with A = and B = represents the charged 
analogue of Tolman's solution VI and thus in the absence of the total charge 
q reduces to the neutral one (the sub-case C of uncharged fluid sphere in Pant & 
Sah^^). The corresponding total gravitational mass m(r = a), can be obtained as 



na 2 (2 -n) + 2q 2 
2(l + 2n-n 2 )a 



(20) 
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4. Detailed study of some of the solutions 
4.1. Bayin's solutions 
Case I: For n — 1 

From the expressions for pressure p(r) and density p(r), as given in Eqs. and 
(|8|). for the central values, viz., r = 0, we have po/po = — 1 which can be written as 
Po + Pa = 0. This equation of state gives the vacuum fluid case. Also this solution 
shows a pressure distribution which goes to negative infinity as r approaches zero 
and density approaches to positive infinity at the centre. According to Bayin^l, this 
solution can be used to represent portions of stars, but can not provide a physically 
acceptable model for the entire star. 




Radius Radius 

Fig. 1. Bayin's solutions showing pressure - density - radius plot for n=l(left) and n=3(right) 
(Case I). 



Case I: For n = 3 

At the origin, r = 0, we get p c /p c = — 1/3 which turns into 3p c + p c = 0. This 
equation of state gives the radiation case. Also the solution for p(r) shows a pressure 
distribution which gives negative infinity pressure at the centre. So it can be used 
also to represent a portion of a star. Here mass density is positive infinity at the 
origin. 

Case II: For n = 1 

At the origin it gives p c j p c = — 1 and we get p c + p c = 0. This equation of state 
represents vacuum fluid case. From the expression of p(r) we see pressure goes to 
negative infinity as r approaches zero. Likewise the above case it can be used to 
represent a portion of the star where mass density is positive infinity at the origin. 
Case II: For n = 3 

At the origin, we get the pressure and density as follows, 



_ 2W + d 



(21) 
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Density Density 

Fig. 2. The pressure vs density (equation of state) plot for n=l(left) and n=3(right) for Bayin's 
solutions (Case I). 



3C 

8 ^ = ( 22 ) 

Here we see that, at the centre, pressure and density have finite values. In order 
to have positive pressure and density at the origin we must have the following 
conditions: C\ < and Wo > |Ci|/2 where C\ and Wo are constants. 
Now, from the Eqs (H]) and (22]) we get 

Pc oGi 6 

This is exactly the same result as Bayin^lgot in his non-charge case. Therefore, it is 
very interesting to note that for both the cases, corresponding to static neutral-fluid 
sphere and static charged-fluid sphere, we get the same result for central pressure 
and density. 

In Bayin's ^31 paper we see p(r) increases as r goes from centre to surface which, 
indeed, is not to be a physically reasonable case. But in our expression we have 
an extra part HK 2 r 2 /8ir which is arising due to inclusion of charge where if is a 
constant. So we can choose suitable values for K such that p(r) can be a decreasing 
function of radius from centre to surface and hence the solution can provide a 
physically valid case. Therefore, for p(r) to be a decreasing function of r from 
centre to surface we must have, from the Eqs. (fT"4|) and (|22|) . the condition 



K - S - I^V) < (-"A ) , (24) 

w 2 wl ) \ wl > 



which gives WqK 2 > 5/11. Thus, altogether we get three conditions as follows: 
Ci < 0, Wo > |Ci |/2 and W§K 2 > 5/11. 

Again, for p(r) to be a decreasing function of radius from centre to surface of 
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the star, we get the condition, from the Eqs. (fTS")) and (f2"Tj) . as 



2 {Cir 2 + W 2 - (1 - 2W 2 K 2 )r i } 
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The above result is obtained on the condition 5/9 < WqK 2 < 1. If we apply this 
condition for radius of the star, for r = a, then we can conclude that through-out 
the star the pressure and density have finite values and they are decreasing function 
of radius from centre to surface. 

Now, we can get the expression for radius of the star considering p(r = a) = 
at the boundary as 
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By matching interior and exterior solutions on the boundary, r = a, we can evaluate 
the integration constants G\ and C2 in terms of mass m and radius a of the star 
which can be expressed as follows: 
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Therefore, after substituting the above value of C\ we can get a relation between 
m and a as 
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The equation of state for this case which is valid throughout the sphere is given by 
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Case III : For n = 1 

When K 2 < 5/9 pressure is negative infinity but density is positive infinity at 
the centre of the sphere. Again, if we take the condition K 2 > 1 then the result 
will be reverse one. We get, from the boundary condition, the value of the constant 
C3 = (1 + K 2 )/m. Equating pressure to zero at the boundary, r = a, we get the 
relation between m and a of the star as follows 



This represents the charge analogue to Bayin's^^case III where the relation between 
the radius and mass was a — m for the value K — 0. Again, for our present charged 
case even putting K — 1 we can recover the same result a — m. In this connection 
it would be very interesting to note that, from the relation q(r) = Kr n , for K = 1 
and n= 1 we get the total charge of the sphere as q — a. As a result, q = m which 
provides the stability condition regarding a charged fluid sphere. On the choice of 
suitable value for K we can get different relation between q and m. A detail study 
of these type of relation a — am have been performed by several worker in different 



In Figs. [T]&E1 we have made a graphical study of Bayin's solutions for Case I, 
setting all the constants to unity, with the only variable being radius, and for the 
different n values. This shows a general trend of the nature of the curve, but fails 
to match with the realistic stars (compact and charged). In a future work, we will 
match conditions to realistic stars, estimating the exact values of the parameters, 
where we can hopefully give a better insight of the solutions. It is also noteworthy 
to mention that not all of the Bayin's solutions are physically feasible to match 
realistic astrophysical objects, nevertheless they are of academic importance. 

4.2. Tolman's solutions 
Case I: For n = 

If we now make the specific choice n — for the parameter n appearing in the 
above solution set then one get the following expressions. 



a = am 
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This gives the equation of state 
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(36) 
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and the mass-radius relation 



(37) 



We see that it gives EMMM for imperfect fluid case. Now if we look at the pressure 
profile, we see it is negative infinity at the centre and increases from centre to 
surface. The pressure has a finite value at the boundary. 
Case II: For n = 0.5 

We are in favor of this case because of some historical reasons. This choice was 
originally done by Tolman himself in his uncharged version. In this case the 
gravitational mass becomes 

3a 2 + 8g 2 

m=^ r . (38) 
The equation of state, by virtue of Eqs. (JT7J) and (|18[) . can be written as 
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which for the present case reduces to 
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In this case also pressure is negative infinity at the centre. But if we put the condition 
(q/a < ±1/3) we see pressure becomes positive throughout the sphere. So, finally 
we get this case as a perfect fluid case having positive pressure but not EMMM. 
At the boundary, equating pressure to zero, we get a 2 = 9q 2 . So, from Eq. (f3"5| one 
can obtain m = (5/6)g which gives q/m > 1. 
Case III: For n = 1 

For this choice of n, the gravitational mass becomes 

a 2 + 2q 2 



4a 



whereas equation of state is 
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(44) 



In this case to have positive pressure everywhere in the sphere we must have the 
condition (q/a < ±1/2). It is also a perfect fluid case having positive pressure but 
not EMMM. At the boundary, equating pressure to zero, we get a 2 = 4q 2 . So, from 
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Eq. (|4"2"|) we have m = (3/4) q which gives q/ra > 1. 

Case IV: For n = 1.5 

In this case the gravitational mass is 

3a 2 + 8g 2 



m 



14a 



and the equation of state becomes 



P + P = 



147TT"' 



2q 2 



The gravitational mass for this case is the same as that of n = 0.5 case. 



1 



567rr 2 
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25g 2 



(45) 



(46) 



(47) 



To have positive pressure throughout the sphere we get the condition (q/a < ±3/5). 
At the boundary, equating pressure to zero, one can get a 2 = (25/9)g 2 . So, from 
Eq. (|4"S"|) we get m — (7/10)g which gives q/m > 1. 




Fig. 3. Tolman's solutions, showing pressure - density - radius plot for n=0.5 (top left), n=1.0 
(top right), n=1.5 (bottom left) & n=2.0 (bottom right) for cases II - V respectively. 
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Density 



O.i 0.12 0.14 



Fig. 4. Equation of states for n=0.5 (top left), n=1.0 (top right), n=1.5 (bottom left) & n=2.0 
(bottom right) for Tolman's solutions for cases II - V respectively. 



Case V: For n = 2 

The gravitational mass in this case becomes 



(48) 



This mass expression is exactly the same as that of the n = case of our previous 
work (Ray & Das, 2004) which vanishes for the vanishing electric charge and thus 
provides 'electromagnetic mass' model (Lorentz, 1904; Feynman et al., 1964). How- 
ever, the equation of state for the present situation differs from that of the n = 
case and is given by 
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2q 2 



(49) 
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To have positive pressure throughout the sphere the condition is q/a < ±2/3. It 
is also an interesting case that this condition gives us it to be a perfect fluid case 
having positive pressure and also it gives EMMM. Otherwise we can generalize all 
above conditions to one condition which is q/a < ±1/3 and by this we get positive 
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pressure for all cases and they are perfect fluid in nature. At the boundary, equating 
pressure to zero, we get q/a — 2/3. So, from Eq. ([48]) we get m = (2/3)q which 
gives q/m > 1. In Figs. [3]&|4l we have plotted the equation of state of the cases 
described in Tolman's solutions for the 4 cases (n=0.5, 1.0, 1.5 & 2) where we have 
set the uniform condition q/a < ±1/3 for all the four cases. For all the above cases 
mass density is positive infinity at the centre and decreases from centre to surface. 
It is again to be stressed that the plots are shown here to give an implication of the 
nature or behavior of the solutions, and matching the variables to real units needs 
some normalization factors, which is not the primary purpose of this paper and will 
be shown elsewhere. 

5. Conclusions 

From discussion in section IV (A) , we do not have any physically viable results for 
Bayin's solutions for case I, n = 1,3 and case II, n = 1. All three cases give negative 
infinity pressure at the centre. In this context we would like to mention that idea of 
negative pressure is not new in astrophysics, especially in the realm of cosmology. 
The equation of state, viz., P + V = for positive mass density provides a repulsive 
type negative pressure MU | j1 | j4jj | j n sense5 negative pressure in our solution 
is not unrealistic. However, for case II, n = 3 we get physically interesting results 
and hence have carrie d o ut detailed study. At origin, the equation of state is the 

|OQ| 

same as that Bayin's with neutral case. Throughout the sphere the equation 
of state is very complicated. Moreover, for Bayin's case pressure was a increasing 
function of radius from centre to surface and hence according to him it was not 
a physically reasonable case. But in our case due to inclusion of charge we get a 
condition by which pressure becomes a decreasing function of radius from centre to 
surface and thus becomes a physically acceptable case. Also, this perfect fluid case 
provides EMMM with positive pressure. Therefore, it gives an example which goes 
in contradiction to Ivanov's 4 - : conclusion that "Electromagnetic mass models are 
subcases, often spoiled by negative pressure" . For this case, we have also found out 
a relation between total mass and radius. For case III, n = 1, we get two conditions: 
one for negative infinity pressure with positive infinity density and the other one 
with reverse features, both at the origin of the sphere. For this case, we get also 
a mass-radius r elat ion which gives again the same result in the absence of charge 
related to Bayin^Slcase. It is interesting to note that for a particular value of K = 1, 
we get total mass of the spherical system equals the total charge. Unfortunately, 
here the fluid pressure becomes negative infinity in the centre. Otherwise, for other 
values of K we get finite values for the ratio of total charge to total mass. 

In the Pant & Sah cases, which are charge analogue to Tolman VI type, 
we have a series of conditions to get positive as well as finite pressure and density 
throughout the sphere. However, from this series of conditions one can easily arrive 
at a general condition which is q/a < ±1/3. One particular case, for n = 2 gives 
EMMM with positive pressure under the perfect fluid condition. Again it is an 
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example which goes in contradiction to Ivanov's conclusion. We also notice that 

for n = 0.5, 1.0, 1.5, 2.0, total charge to total mass ratio becomes constant in each 

case and it is greater than unity. It reminds the general theorem given by de Felice 
1041 

et al. that if the total electric charge of a perfect fluid ball is smaller than its 
total mass, then there is no regular static configuration having a radius arbitrarily 
close to the size of the external horizon. Therefore, the above four cases again give 
support in favor of the stability of the charged spherical models in connection to 
normal stars. 
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